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The Tensor- Vector-Scalar theory of gravity, which was designed as a relativistic implementation 
to the modified dynamics paradigm, has fared quite well as an alternative to dark matter, on both 
galactic and cosmological scales. However, its performance in the solar system, as embodied in the 
post-Newtonian formalism, has not yet been fully investigated. Tamaki has recently attempted to 
calculate the preferred frame parameters for TeVeS, but ignored the cosmological value of the scalar 
field, thus concluding that the Newtonian potential must be static in order to be consistent with the 
vector equation. We show that when the cosmological value of the scalar field is taken into account, 
there is no constraint on the Newtonian potential; however, the cosmological value of the scalar 
field is tightly linked to the vector field coupling constant K, preventing the former from evolving as 
predicted by its equation of motion. We then proceed to investigate the post-Newtonian limit of a 
generalized version of TeVeS, with jEther type vector action, and show that its /3,'y and ^ parameters 
are as in GR, while solar system constraints on the preferred frame parameters ai and 02 can be 
satisfied within a modest range of small values of the scalar and vector fields coupling parameters, 
and for values of the cosmological scalar field consistent with evolution within the framework of 
existing models. 

PACS numbers: 04.25.Nx, 04.50.Kd, 04.80.Cc 



I. INTRODUCTION 

As is known. General Relativity (GR) cannot explain 
the dynamics of our universe on large physical scales, 
since the amount of visible mass clearly lies below what 
would be expected from the observed gravitational ef- 
fects. The usual remedy is to invoke a form of matter 
which does not couple to light, therefore being referred 
to as Dark Matter (DM). However, one can also take a 
different point of view and modify the law of gravity it- 
self. Such a solution to the missing mass problem was 
first studied in great detail by Milgrom in his MOND 
paradigm. 

The modified Newtonian dynamics (MOND) 
paradigm [T], proposes that Newtonian gravity progres- 
sively fails as accelerations drop below a characteristic 
scale ao — 10~^°m/s^ which is typical of galaxy out- 
skirts. MOND assumes that for accelerations of order Oq 
or well below it, the Newtonian relation a = —W^n is 
replaced by 

/i(|a|/ao)a= -V$^v, (1) 

where the function /i(a;) smoothly interpolates between 
//(x) — X ai X 1 and the Newtonian expectation 
jl{x) = 1 at a: 1. This relation with a suitable stan- 
dard choice of jl{x) in the intermediate range has proved 
successful not only in justifying the asymptotical flatness 
of galaxy rotation curves where acceleration scales are 
much below ao, but also in explaining detailed shapes of 
rotation curves in the inner parts in terms of the directly 
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seen mass, and in giving a precise account of the observed 
TuUy-Fisher law which correlates luminosity of a disk 
galaxy with its asymptotic rotational velocity This 
sharp relation, while obtained naturally in the framework 
of MOND, requires quite a fine tuning of dark halo pa- 
rameters to be explained by the dark matter paradigm. 

However, MOND alone is only a phcnomenological pre- 
scription that does not fulfill the usual conservation laws, 
nor does it make clear if the departure from Newtonian 
physics is in the gravity or in the inertia side of the equa- 
tion F = ma. Moreover, it is non relativistic, and as 
such it does not teach us how to handle gravitational 
Icnsing or cosmology in the weak acceleration regimes. 
To address these issues, Bekenstein designed TeVeS [S], 
a covariant field theory of gravity which has MOND as 
its low velocity, weak acceleration limit, while its non- 
relativistic strong acceleration limit is Newtonian and its 
relativistic limit is general relativity (GR). TeVeS sports 
two metrics, the "physical" metric on which all matter 
fields propagate, and the Einstein metric which interacts 
with the additional fields in the theory: a timelike dy- 
namical vector field. A, and a scalar field, (p. The theory 
also involves a free function J-', a length scale £, and two 
positive dimensionless constants k and K. The scalar 
field in TeVeS provides the additional gravitational po- 
tential for matter, whereas the vector field provides the 
desired light bending properties, in a fashion similar to 
the constant unit vector in Sanders' stratified theory [4]. 

Many aspects of TeVeS have been investigated exten- 
sively, proving the theory to be faring quite well in view 
of the huge challenges it was designed to meet. Beken- 
stein showed that TeVeS's weak acceleration limit repro- 
duces MOND, and that it also has a Newtonian limit [3]. 
Skordis [S| formulated the cosmological equations for 
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TeVeS, for both background and linear perturbations, 
and later extended his investigation to a version of TeVeS 
where the action of the vector field is of Einstein- iEther 
form [B^. The CMB spectrum and the matter power 
spectrum P{k) were calculated by Skordis, Mota, Fer- 
reira and Boehm [:7|, who showed how TeVeS can repro- 
duce the power spectrum in a manner similar to Dark 
Matter. Further inquiries into TeVeS cosmology have 
been made by Dodelson and Liguori [5] , who showed that 
perturbations in the TeVeS vector field can drive struc- 
ture growth, and by Bourliot et al., [3] who considered 
a broad family of functions that lead to modified gravity 
and calculated the evolution of the field variables both 
numerically and analytically. Giannios |10j has found 
exact solutions of TeVeS for spherically symmetric sys- 
tems, including Schwarzschild-like black holes, and Sagi 
and Bekenstein expanded upon his work and found 
charged black hole solutions. Zhao and Famaey [12] have 
put a variety of constraints on the TeVeS free function 
from galaxy dynamics. Laski, Sotani and Giannios [T3] 
investigated neutron stars in TeVeS, using them to place 
a lower bound on the allowed cosmological value of the 
scalar field, and Sotani [2] calculated the fundamental 
oscillation modes of neutron stars for the theory, showing 
how the imprint of the scalar field could be detected in 
gravitational waves. 

TeVeS has also been tested against a multitude of data 
on gravitational lensing. Chiu, Ko and Tian have exam- 
ined theoretical predictions of TeVeS for amplifications 
and time delays in strong gravitational lensing [15j , while 
Zhao et al. have put TeVeS predictions for image spHt- 
tings and amplifications to test against a large sample of 
lensed quasars [16_. And Chen and Zhao have compared 
the statistics of strong gravitational lensing by galaxies 
with TeVeS [T7| . This work is admirably capped by Chen, 
who calculated the lensing probability with image sepa- 
ration larger than a given value Ad in an open, TeVeS 
cosmology, and showed that the predicted lensing prob- 
abilities with the simple interpolating function x/{l + x) 
match the observational data quite well [TH]. Angus et 
al. have criticized the claim that the colliding clusters of 
galaxies "the bullet" pose a threat to gravitational lens- 
ing a la TeVeS [IS]. 

However, it is not yet clear where the theory stands 
with respect to solar system constraints, which are usu- 
ally embodied in tight limits on the allowed values of 
post-Newtonian (PN) parameters. Any general metric 
theory of gravity can be fully characterized by ten pa- 
rameterized post-Newtonian (PPN) parameters PTT, that 
quantify the lowest order effects in w^/c^ and GnM /c^r. 
Five of these parameters, C11C2, Csi C4i ^^'^ Q^3i vanish 
identically for any semi-conservative theory, i.e. one de- 
rived, like TeVeS , from a covariant action principle. Two 
others, known as the Eddington— Robertson— Schiff pa- 
rameters (3 and 7, characterize respectively the nonlinear- 
ity and the spatial curvature produced by gravity. Of the 
remaining three PPN parameters, two, ai and a2 charac- 
terize preferred frame effects, and the third, ^, also known 



as the Whitehead parameter, characterizes a peculiar sort 
of three-body interaction. TeVeS' PPN parameters have 
been calculated only under simplifying assumptions, such 
as spherical symmetry [3J [TD] or cosmological value of 
the scalar field set to zero [3T]. These assumptions did 
not enable calculation of the preferred frame parameters 
ui and a2- In a theory possessing a preferred Lorentz 
frame, the rest frame of the timelike vector field, one ex- 
pects these parameters to be different from zero, whereas 
solar system experiments |22j constrain their measured 
value to be very close to zero, specifically ai < 10^^ 
and 1 02! ^ 4 x 10""^. The question then rises, is there a 
region of TeVeS parameter space within which the pre- 
ferred frame parameters are zero as in GR, or is TeVeS 
ruled out or unreasonably constrained by solar system 
tests? 

In this work, we start with an overview of TeVeS in 
Section |TT] and elaborate on the PPN parameters in Sec- 
tion |III| giving the full set of PPN parameters for TeVeS 
with no prior assumptions such as spherical symmetry 
or zero cosmological value of the scalar field. Although 
we find (3 and 7 to be unity, and ^ = 0, as in GR, we 
encounter a difficulty with the vector equation, and find 
that to PN order, it links the cosmological value of the 
scalar field to the coupling parameter of the vector field K 
in such a way that the scalar field is not allowed to evolve 
with cosmological expansion. This is another hint to the 
fact that the simple Maxwell- like action of the vector field 
causes dynamical problems: it has previously been shown 
by Seifert i23j that static spherically symmetric solutions 
in this theory are unstable against spherically symmetric 
perturbations, and by Contaldi et al. \2S that the vector 
field runs into caustic singularities in rather generic situ- 
ations. Therefore in Section [TV] we present results for the 
PPN parameters for a more general form of TeVeS , one 
with a full vector action of iEther type. We show that the 
preferred frame parameters can be set to zero within a 
modest range of small coupling parameters of the scalar 
and vector fields, and for reasonable values of the cos- 
mological scalar field. We conclude that for particular 
ranges of the coupling parameters, TeVeS with general- 
ized vector action is indiscernible from GR in the solar 
system. The full details of the calculations are given in 
the appendices. 



II. TEVES 

TeVeS has MOND as its weak potential, low acceler- 
ation limit, while its weak potential, high acceleration 
limit is the usual Newtonian gravity. TeVeS is endowed 
with three dynamical gravitational fields: a scalar field 
(j), a timelike unit normalized vector field A", and the 
Einstein metric gap on which the gravitational fields of 
the theory propagate. The theory also employs a "phys- 
ical" metric (jaiB on which gauge, spinor and Higgs fields 
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propagate. It is related to ga[3 by 



9a/3 = e '^'^gaf3 - 2AaApsmh.{2(l)). 



(2) 



The index of or of </> ^ is always raised with the metric 
g"^ , the inverse of gap. 

The equations of motion for the fields in TeVeS derive 
from a five-term action depending on four parameters: 
the fundamental gravity constant G, two dimensionless 
parameters k and K and a fixed length scale i. The 
familiar Einstein-Hilbert action for the metric and the 
matter action for field variables collectively denoted / 
have the form 



-gd X, 



Sr, 



(3) 
(4) 



Next comes the vector field's action, with K a dimen- 
sionless positive coupling constant 



K 

'32ttG 



-^{g^'^A.A^ + l) 



(5) 



which includes a constraint that forces the vector field to 
be timelike (and unit normalized); A is the corresponding 
Lagrange multiplier. The presence of a nonzero A" estab- 
lishes a preferred Lorentz frame, thus breaking Lorentz 
symmetry in the gravitational sector. Finally, we have 
the scalar's action (k is a dimensionless positive param- 
eter while £ is a constant with the dimensions of length, 
and a dimensionless free function) 
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T{kt^h'^^(j,,a<j>,p) V^d^ 



(6) 



Above /i"'' = g^f^ - A°'Af^ with A" = g^f^Ap. 

Variation of the action with respect to g"'^ yields the 
TeVeS Einstein equations for g^p 



Gal: 



Tafj 



(7) 

where t'(Q,A^) = v^Ap + A^vp, etc. The sources here are 
the usual matter energy-momentum tensor T^p (related 
to the variational derivative of S,n with respect to g"^), 
as well as the energy-momentum tensors for the scalar 
and vector fields. 



2k'^PG 



^ap = K (^^''' A^^,^a]A[y^p]~ -g"^ g^"" A\^„^^,}A^r,^]9aP^ 

- ^AaAp 
where V[aAp-\ = VaAp — AaVp, etc., and 



(9) 
(10) 



Each choice of the function T{y) defines a separate TeVeS 
theory. Its derivative ^(y) functions somewhat like the 
fl function in MOND. For y > 0, ii{y) ~ 1 corresponds 
to the high acceleration, i.e., Newtonian, limit, while the 
limit < n{y) ^ 1 corresponds to the deep MOND 
regime. We shall only consider functions such that > 
and /i > for either positive or negative arguments. 

The equations of motion for the vector and scalar fields 
are obtained by varying the action with respect to 4> and 
Aa, respectively. We have 

[/i(y)/i"^0, a] . ^ = fcG [g^P + (1 + A^Af>] f^p , 

(11) 

for the scalar and 



Stt 



= SttG (1 - e-"*-^) g'^'^Af^f^p 



p^XA" + ^^^AU^pg'''^(|)^^ 



(12) 



for the vector. Additionally, there is the normalization 
condition on the vector field 



A"A^ 



9ap 



(13) 



The A in Eq. (12 1, the lagrange multiplier charged with 



the enforcement of the normalization condition, can be 
calculated from the vector equation. 

The three parameters, fc, K and ^, all specific to TeVeS, 
are constant in the framework of the theory, as is G, the 
fundamental gravitational coupling constant, which does 
not coincide with Newton's Gjy. 



III. THE PPN PARAMETERS FOR TEVES 

In the weak field, slow motion limit, the next-to- 
Newtonian order gravitational effects of any metric grav- 
itational theory can be described in terms of a set of 
functionals of the matter variables that answer certain 
criteria of "reasonableness" and simplicity (see [5D] for 
details), known as the Post-Newtonian potentials, and 
of ten parameters, 7, (3, Ci, C2, Ca, C4, oil, "2, as, 
known as the PPN parameters. Since TeVeS has two 
metrics, the Einstein metric on which the gravitational 
fields propagate, and the physical matter metric, serving 
as background for massive and massless matter particles, 
the potentials are of course to be expressed in terms of 
physical coordinates defined by the physical metric, and 
physical fluid variables, with indices raised and lowered 
by the physical metric. 

The standard form of the Post Newtonian metric is 
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given by 
500 = 



-1 + 2U- - 2£,<i>w+ 

+ (27 + 2 + a3 + Ci-20«'i 
+ 2(37 - 2/? + 1 + C2 + 0*2 
+ 2(37 + 3C4- 20*4 -(Ci- 



where is 0(1) and z;* are 0(5). This form of the fluid 
velocity is normalized with respect to the physical metric 



a 13 



-1 



(26) 



ffOj 



^(47 + 3 + al 
1 



(1 



0^2 + Ci 



+ 2(1 + C3)'i>3 

20-4 



gjk ={l + 2-iU)5jk 
The potentials are all of the form 



il(a;, t) = 



_ f p{x',t)f{x',t) 



;3 / 



(14) 

(15) 
(16) 

(17) 



Following the procedure described in ^20', one can 
calculate the post-Newtonian physical metric for TeVeS. 
Since the calculation is straightforward but tedious, we 
relegate it to the appendix. We only mention here, that 
with nonzero cosmological value of the scalar field, the 
physical metric is not asymptotically Minkowski, and 
to be able to compare our result to the standard form, 
we must transform the physical metric to local quasi- 
Cartesian coordinates, through the following coordinate 
transformation 



(27) 



where f{x',t) is given for each potential as follows 
<i>2 ■■ U, <i>3 : n, <i>4 



U -.1, $1 : v'v\ 
$w : / d^x"p{x",t)' 



[vi{x - x')i 



r'|2 



{x-x"), 



P 
P 

{x' - x")^ _ 

\x — x"\ \x' — x'' 
Vj{x — x')j{x — x'Y 



The potentials must then be rescaled accordingly, for ex- 
ample U = e-^'^oU, % = e-'^'^^Vi and Wi = e'^'^'^Wi. 
Additionally, we will work in units such that Gjv = 1. In 
(18)these units, the coupling constant G can be expressed in 
terms of k [25] 



G 



1 



l-K/2 



k 

in 



(28) 



The standard post Newtonian gauge is such that all 
dependence on x,oo and x.ij ti^s been eliminated from 
goo and g^, where x is the "superpotential" defined as 



In local quasi-Cartesian coordinates, and with Gn 
the post-Newtonian metric for TeVeS is given by 



ffoo 



-l + 2[/-2C/2+4$i+4$2 + 2$3 + 6$4, (29) 



X = - / pix',t)\x-x'\d''x' 



(19) 



Thus the coordinate frame is determined up to the 
necessary order, and the form of the metric components 
is unique. 

When taking into account the cosmological value of the 
scalar field, 4>o, the Einstein metric and fields to post- 
Newtonian order have the following form 



gij =Sij{l + 2U) 



(30) 



and 



9o^ = -^{7 + ai ~ a2)V^ ~ ^{1 + a2)W^ (31) 



with 



ai 



4 



,(1) 



(2) 



goo — "-00 

9tj = {l + hiJ)6^j 
goi = hoi 

A°' ^{1 + + A^'^^\A'',Ay, A^) 
(1)^(1)0 + ^^'^ + ^^^^ 



(20) 
(21) 
(22) 

(23) 
(24) 



4G 

Ik 



{{2K ~ 1 



(32) 



and 



Q!2 



2G 



{2-Ky 



(3(2-i^)-(X + 4)e'*'^«)-l (33) 



where an index of (1) denotes a term of first order in 
the dimensionless gravitational potential GnM/c^t (and 
thus second order in the velocity v/c), an index (2) de- 
notes a term of second order in the potential, and A^ is 
of 0(1.5) in the potential. 0o is to be determined by the 
cosmological boundary conditions. The fluid velocity is 
given by 



-'^''{l + v\v'',v\v''), 



(25) 



Here we omitted the bars over the potentials. The 
time-time and space-space components of the physical 
metric do not depend on the cosmological value of the 
scalar field in the standard post-Newtonian coordinate 
system, yielding /3 = 7 = 1, ^ = Ci = a3 = 0, with 
i= 1..4, as for GR. 

The result of /3 = 1 is apparently inconsistent with 
Giannios' [10], who obtained /3 ^ 1 for A'' ^ 0. However, 
Giannios performed the calculation assuming spherical 
symmetry, and took the radial component of the vector 
field to be of 0(1), and not 0(1.5) as dictated by the 
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Post-Newtonian formalism. While such a term is allowed 
by the field equations, the vector equation to 0(1) is 



(34) 



This means that such a term in the vector field does not 
originate in sources or fields, and hence it is not clear how 
it would be incorporated in the PPN formalism. This was 
not noticed by Giannios since he performed the calcula- 
tion in vacuum, so that the sources did not appear di- 
rectly in the equations, but only as boundary conditions 
for integration. An 0(1) term in the vector field would 
require extra parameters to set the boundary conditions, 
and thus seems unnatural. An which is 0(0.5) could 
be acceptable, but it would have no effect on the calcula- 
tion of the PPN parameters, since it could be eliminated 
by an appropriate Lorentz transformation. 

However, the time-space component of the metric does 
depend on ^o, as do the preferred frame parameters, and 
factors of e^'^° cannot be "rescaled out" by a coordinate 
transformation. Therefore it is incorrect to perform the 
PPN parameters calculation without taking into account 
the cosmological value of the scalar field. In fact, its role 
in the vector equation is quite critical, as we shall soon 
see. 

The covariant divergence of the vector equation, unlike 
the covariant divergence of the Einstein equation, is not 
automatically satisfied, but yields a constraint on the di- 
vergence of the vector field or on the coupling constants 
of the theory (see [20] j Section 5.4 for details). In the 
case of TeVeS , we obtain to 0(1.5) 



K 



U,o = -2(1 



(35) 



Here Tamaki's calculation [21] was in error; since he as- 
sumed (po = 0, his conclusion at this step should have 
been that necessarily K = 0, namely one cannot take 
into account the coupling of the vector field in this par- 
ticular theory without the scalar field, because obviously 
one cannot demand C/.o = always. When a nonzero cos- 
mological value of the scalar field is taken into account, it 
is constrained by the divergence of the vector equation. 
From the requirement that U,q ^ 0, we obtain 



K 



1 -K/2 



-2(1 



-400 



(36) 



which links the cosmological value of the scalar field to 
the coupling constant of the vector field as follows 



(37) 



This contradicts the basic assumption of TeVeS that 
its coupling parameters are constant, namely that they 
do not change during the evolution of the universe. Or 
vice- versa, if the assumption remains valid, then 0o is not 
allowed to evolve. Moreover, this relation would mean 
that (/)o is negative, since it has been shown [11 that 
one must have K < 2. This would bring the problem of 
superluminal propagation of scalar waves. 



When relation (37) holds, it can be shown that there is 



no ambiguity in the determination of the physical metric 
and the preferred frame parameters, although the spa- 
tial divergence of the vector field remains indeterminate. 
However, it looks like the constraint (37 1 is another sign 



of the existence of a dynamical problem with the vector 
action in TeVeS , adding up to the analysis of Contaldi 
et al. [21] regarding the formation of caustic singularities 
in the evolution of the vector field. Therefore, we leave 
aside simple TeVeS , and proceed instead to calculate the 
PPN parameters of TeVeS with Jjlther [26] type vector 
action. 



IV. TEVES WITH jETHER TYPE VECTOR 
ACTION 



This version of TeVeS has been proposed by Skordis [6] , 
who investigated its cosmology, as a natural generaliza- 
tion, and Contaldi et al. ,24] adopted it as a resolution 
to the problem of vector caustic singularities. The vector 
action is taken to be of the most general form quadratic 
in derivatives of the vector fields, whereas its scalar and 
metric actions are unaltered, thus preserving the correct 
MOND and Newtonian limits. We take a vector action 
of the form 



5„ = I V^d^x [^F^pF^P + ^S^pS""" + K2 {VAf + KiAo^A" - A (A"A„ + 1)^ (38) 



r 



where Fap = Aa-p - Ap.a, Sap = Aa-p + Ap-a and 2K+,C2 = K2 and C4 = -^4. 

A"^ = A^A'^p. The relation between the coupling con- The equation of motion and the stress tensor for the 

stants Ki and the q of iEther [26] is C1—C3 — 2K,ci+cs = vector field are then 
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= SnG (1 - /"f^^A^ 

0^/5 = A' - \^^9ap) + ^+ {s^aSp" ~ Qo^p + V. [A'' S^f3 - S%^A 



(39) 



A^A^Ap 



9af3 



A^A" - XA^A 



(40) 



and the scalar equation remains unaltered. To post- 
Newtonian order, the Einstein and physical metric and 
fields are as in Section |III| and the coordinate rescaling 
is as for simple TeVeS. The new gravitational coupling 
constant Gn is given by 



The PPN metric is 



500 = -1 + 2C/ - + 4$i + 4$2 + 2$3 + 6$4, (43) 



Gat = G I ( 1 



k 

4tt 



(41) 



g^J ^ S,j{l + 2U) (44) 



and for units in which Gat = 1, 



G= 1- 



K + K+- Ki 



k 

47r 



(-1) 



Att{2 - {K + K+) + Kj) 
' St: + k{2 - {K + K+) + Ki) 



(42) 



and 



90, = -^{7 + ai~ a2)V^ - ^(1 + "2)^, (45) 



with 



ai 



a2 



f G ((e^-^-o - (1 - 2i4:)e-2'^") sinh (2(/)o) - (i^+e*'^" + K)) \ 
^ \ 2KK+ -{K + K+) "^J 

oi 4G ((2 + 37^2 + 2K+)e^^" - {2 - {K + K+ - K^)) er^^") sinh (20o) 

{K2 + 2K+) {2~{K + K+- Ki)) 
2G ((if4 -K + iK+){2 + 3X2 + 2is:+)e'^'^« + (3ii:+ + 3/^2 + X - Ki){2 - [K + 



-Ka))) 



{K2 + 2K+) {2~{K + K+- Ki)Y 



(46) 

(47) 
(48) 



Hence TeVeS with iEther type vector kinetic term has 
/3 = = 1, ^ = = as = 0, with i — 1..4, as for GR, 
and preferred frame parameters ai and a2 given above. 

Although not apparent from the above formulation of 
the alphas, when the coupling constants of TeVeS, k and 
the Ki, as well as the cosmological value of the scalar 
field are zero, then ai — a2 = Q- Moreover, when the 
scalar field is decoupled from the theory by setting k = 
(j)o — 0, the preferred frame parameters acquire the values 
calculated for jEther theory, as given in [27] 



ai 



2KK+ -{K + K+) 



(49) 



ai (K - 3K+ - Ki)(K + SK+ + 3K2 - KA 

(Xo — h 

2 {2K+ + K2){2 - {K + K+ - Ki)) 

(50) 

This was to be expected, since when TeVeS is exempted 
of the scalar field, the physical and Einstein metric coin- 
cide, and the theory is identical to jEther. 

Current data strongly constrains the preferred frame 
parameters to |ai| < 10^* (from measurement of the 
earth-moon orbital polarization with lunar laser ranging 
and from pulse timing of binary pulsar PSR J2317-I-1439) 
and \a2\ < 4 x 10~^ (from solar alignment with the eclip- 
tic) gHllSllin]- It has been shown that in jEther theory 
the two alphas can be set to zero with two parameters 
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to spare ^Tj; it is obvious that in the extended version 
of TeVeS, which has five parameters, we can solve for 
ai,a2 ~ with three parameters to spare. However, 
since the alphas depend on the cosmological value of the 
scalar field, their value is expected to change with cosmo- 
logical evolution, and the question that rises is whether 
the present day cosmological value of the scalar field, as 
evolved from reasonable initial conditions, is consistent 
with the experimentally measured values of ai and a2 
today, and with small values of the coupling constants fc, 
Ki. Small values are required since in simple TcVcS , in 
order to avoid too large MOND effects in the solar sys- 
tem, k should be of the order of 0.01, and in order to give 
correct lensing K should also be small and of the same 
order of magnitude as fc [T. It then seems natural to 
assume that the extra kinetic terms added to the vector 
action should also have small couplings. 

There is also the question of cosmological evolution 
of the alphas themselves; however, since all experimental 
data on the values of the preferred frame parameters orig- 
inates nearby by cosmological standards (solar system or 
pulsars in our galaxy) , we allow ourselves to assume that 
4>o is constant for practical purposes. 

A. Allowed Ranges of Parameters 

Solving ai ^ and a2 = 0, we obtained expressions 
for K2 and k in terms of K4 , K, and 0o • The con- 



Both k and K2 no longer depend on 7^4, and K2 doesn't 
depend on (pQ. However, Skordis [6] has shown that cos- 
mologically the combination K~\-K+ —K4 in TeVeS with 
iEther type vector action plays the same role as K in 
simple TeVeS, by allowing for a growing mode in the 
vector field in order to source structure formation. For 
this same reason we have not considered negative values 
of K + K+ - Ki. 



V. CONCLUSIONS 

We calculated the preferred frame parameters for 
TeVeS, and in the process of calculation we found that 
the divergence of the vector equation constrains the cos- 
mological value of the scalar field to be related to the 



straints on the alphas are satisfied for large ranges of the 
parameters K^., K, K^; however, when demanding that 
the parameters be small with respect to unity, the ranges 
narrow significantly. One can look at the allowed region 
in the K — A% parameter space, obtained from the over- 
lap of the regions in which both k and K2 are small, for 
specific values of the parameter K4. We take 0o = 0.003, 
as suggested by jl3j, and show, for example, the region 
in the K — plane for which ai,a2 — 0, K2 and k 
are in the range [0,0.3], with K + — K4 = c, for 
c = 0.01 and c = 0.001, in Fig.Q and Fig.(|t|, respec- 
tively. The allowed ranges for small parameters are 
quite narrow, especially for K+, which is restricted to a 
small range near 0.01. The plots remain similar if in- 
stead of requiring ai = and a2 = 0, we only demand, 
for instance, ai = 10^^ and a2 = 0, in compliance with 
current experimental constraints. 

If we increase (j)Q, the shape of the allowed region in 
the K — plane remains similar, but it is shifted up- 
wards, favoring larger values of K^. For too large a value 
of 00 (around 0.01), there would be no overlap between 
the regions in which both k and K2 are small; the same 
happens when K + — K4 is taken to be too large. 
When K + — K4 = 0.1 there is no overlap between 
the two regions. 

In the special case K + Kj^ — K4 = 0, the forms of if 2 
and k become particularly simple 



(51) 



(52) 
I 

coupling parameter of the vector field. Since such a link 
does not allow the scalar field to evolve cosmologically, 
we set aside simple TeVeS and proceeded to calculate 
the preferred frame parameters for a generalized version 
of TeVeS previously suggested in the literature to resolve 
possible dynamical problems [6l El] , TeVeS with jEther 
like vector action. We obtained expressions for ai and a2 
in terms of the couphng parameters of TeVeS , k and Ki, 
and the cosmological scalar field (j)o . Since all existing ex- 
perimental data on values of the alphas originates from 
within our galaxy, we assumed the cosmological value 
of the scalar field to be constant, and of modest value, 
and analyzed the allowed ranges of coupling parameters 
for which both preferred frame parameters are zero. We 
found that the conditions ai = 0, 012 = can be satisfied 
for small ranges of the coupling parameters, and in these 



2tt {{2K - l)e-^^" - {2K+ - l)e*'''» - 2{2KK+ - {K + K+) + 1)) 
2KK+ -{K + K+) 

^ -2{K+{2K-l)+2K) 
^ 3(2if-l) 
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FIG. 1: The region in the K — plane for which ai, a2 ~ 0, 
with 00 ~ 0.003, k and K2 in the range [0,0.3], and K + 
K+ — K4 = 0.01 is shown in color. The shading of the region 
shows the value of K2; a lighter hue corresponds to a higher 
value. K2 increases with increasing K, whereas k increases 
with increasing 




FIG. 2: The region in the K — Kj^ plane for which ai, 122 = 0, 
with 00 = 0.003, k and K2 in the range [0,0.3], and K + K+ — 
Ki = 0.001 is shown in color. The shading of the region 
shows the value of K2; a lighter hue corresponds to a higher 
value. K2 increases with increasing K, whereas k increases 
with increasing K+. 



ranges TeVeS has its PPN metric identical to GR, mak- 
ing it indiscernible from the latter in the solar system. 
Future work on gravitational wave speeds and stability 
is expected to constrain the allowed values of parameters 
further. 
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APPENDIX A: CALCULATION OF PPN 
PARAMETERS FOR TEVES 

The PPN parameters are determined by solving the 
gravitational field equations with a perfect fluid source, 
in a standard coordinate gauge where the spatial part of 
the metric is diagonal and isotropic. The fluid variables 
are assigned orders oi U ^ v'^ ^ p ^ H ^ p/ p ^ 0{1), 
where U is the Newtonian gravitational potential. Taking 
the time derivative of a quantity raises its order by one 
half, since d/dt ~ Vw. In the Newtonian limit of any 
theory of gravity, the temporal component of the metric 
500 is required to first order in the potential, whereas the 
spatial components of the metric gij are of zero order; the 
post-Newtonian limit calls for knowledge of t^oo to 0(2), 
gij to 0(1), and 170^, which must change sign under time 
reversal, is order 0(1.5). The standard form of the Post- 
Newtonian metric is given in Eqs.(14)-(16l. 



In the calculation we will use the following relation 
that holds for U, <i>i,2,3.4, and Vi 



V^ii=~iTrGpf. 



The superpotential Eq.(19l satisfies 

-2U 



-Uo 



(Al) 



(A2) 



(A3) 



(A4) 



the last two relations follow from the continuity equation 
for the fluid, assumed to hold to (9(1.5) 

P.o + {pv')a = 0. 

The standard post Newtonian gauge is such that all 
dependence on x,oo and x.ij has been eliminated from 
500 and g.y . 

We will follow the procedure described in Will's clas- 
sic reference [20] • Although many details are similar to 
Tamaki's computation |21j . we will show the calculation 
again to make it clear where factors of e^^"^" enter, and 
how they alter the results. In TeVeS, which is endowed 
with a scalar and vector field additional to the metric, 
they must be solved for to the order necessary to derive 
the metric as required; namely, 4> to 0(2) in the Newto- 
nian gravitational potential, A* to 0(2) in the Newtonian 
gravitational potential, and to 0(1.5). Moreover, the 
metric from which the PPN parameters are to be deter- 
mined is the physical metric, gap, since this is the metric 
on which matter fields propagate. For the calculation it 
is convenient to rearrange the Einstein equations ([?]) as 
follows 



R. 



inG (f^„ + (1 - 6-4"^) A-^f^c^A,) + T,,,) 



9,.] - i^g^pgn 



(A5) 
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where Tap and 0^/3 are given by Eqs. ([8|-([9| respectively. 
For a perfect fluid source, we take as customary 



Tap = {p+ pTl+p) UaUi3 + pgap 



(A6) 



In the Solar system, where accelerations are high com- 
pared to the MOND acceleration scale gq, we can take 
/i « 1, and, correspondingly, J-{y) ~ y. 



With the Einstein metric and fields of the form (20) 



( 24 1 , the physical metric to the required order is given 

by 



500 = 6^^" (-1 + 
'4o^-20(^) 



(4o^ - 20(^ 
-2(0«)2 + 2<^(i)/iW 



00 



9ij 



= ( 1 + ( h,, - 2(i(i) 



go, = e^'^o/ioi + 2 A' sinh 20o 



(A7) 
(A8) 
(A9) 



where different orders were separated by brackets. 

In our calculation, we will impose the following gauge 
conditions 



) 

00 J- 



hk 



kk,j 



hok.k 



2-K 



Ikk.O 



(AlO) 



(All) 



the first condition is frequently used in GR, whereas the 
second condition is required to bring goo to the standard 
post-Newtonian gauge. Here and throughout the calcu- 
lation we assume the Einstein summation convention to 
relevant order, meaning that repeated spatial indices are 
summed over. Indices on the vector field are raised and 
lowered with the Einstein metric, and indices for the mat- 
ter fields and fiuid velocity are raised and lowered with 
the physical metric. 



1. Calculation of A' and from normalization 

We are now ready to proceed with the calculation. 
First, we find the temporal components of the vector field 
and of the fiuid velocity from the normalization condi- 
tions Eqs. (13 1, (26), respectively 



^ — 2 00 



.t(2) _ 1,(2) 3 (1) 2 

^ — 2 00 8 



2 



20 



(1) 



(A12) 



(A13) 



(A14) 



2. 500 to 0(1) 
Next, we solve for goo to 0(1). To this end, we must 



find h'^od and d)'-^^ 



stein equations (A5) to 0(1) is 
1 



The temporal component of the Ein- 

(A15) 



(1 - K/2) V^^o = -47rGe^20o 



yielding for hf^J 



'oo 



2Ge-2'^'' 



P 



1 - K/2 J |x - x' 
For convenience, we define 

G 



P 



(A16) 



G 



K 



1 - K/2 



(A17) 



since this is a recurring quanti ty i n the calculation. 
6^ ' is determined from Eq.(hl| to 0(1) 



givmg 



kGe-'^'t'" 
47r 



P 



(A18) 



(A19) 



Combining Eqs. (|B2p,(|A19| and (|A7l we have 

P 



goo = Ge 



-200 



1 



1 - K/2 4tt J J |x - X 



-d-^x' (A20) 



We now put the physical metric in standard Newtonian 
and post-Newtonian form in local quasi-Cartesian coor- 
dinates, as described in Ref. [20J, through the following 
coordinate transformation 



e'l^'x^ ,x' = e-^^x' 



(A21) 



In the new coordinates, the physical metric to Newtonian 
order is 



ffoo = -1 + 2GnU , (jij = Sij ,goi = 



with 



Gjv — G 



1 



k 



1 - K/2 47r 



(A22) 



(A23) 



Here the potential U has been redefined in accordance 
with the coordinate transformation U — e~'^'^°U. In the 
asymptotically Minkowskian coordinates, we have 



h'i^ = 2GkU 



and 



kG 

4-77 



u 



(A24) 



(A25) 
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In order to compare our result with ( 14 1 , we must select 



units such that Gn = 1. The coupling constant G is then 
given in terms of K, k 



G 



1 



(A26) 



(A27) 



I - K/2 47r^ 
so that finally 

5oo = -1 + 2C/ , gij = 6ij ,goi^O 

3. g,j to 0(1) 



The next step will be to solve for gij to 0{1). From 
the space-space component of Eq. (|A5[) we have for hij 



hij — Hqq. ij +hkk,ij — hki J k~ hkj ,ik 



making use of the gauge Eq. (AlOl, this becomes 



-8^TGKe~^''"'pS^j, 
(A28) 



whose solution is 

hij — Sij2GK 



or 



2GA'e 



-200 



U6, 



(A29) 



(A30) 



(A31) 



Combining this result with Eq. (A19l we obtain for 



g,j = e"2^"%(l + 2e-2^«Gjv;7) 



(A32) 



We now transform the coordinate system according to 
( A21 1, set units such that Gjv = 1, and obtain 



g,, = <5,j (1 + 2U) . 



(A33) 



Comparing this with Eq. ( |16| we see that for TeVeS 
the PPN parameter 7=1. This is in accordance with 
previous work [3J (TUl [51] , and agrees with current solar 
system data. 



4. go, to 0(1.5) 

We now solve for goj to 0(1.5), which requires calcu- 
lating and h^j to the same order. First, we find the 
Lagrange multiplier A to 0(1) from the temporal compo- 
nent of the vector equation Eq. ( 12 ) 



A -A'V^^o - 167rGpsinh20o (A34) 



The time-space component of Eq. (A5l to 0(1.5) is 
then 



iQk.kj + hkk.Qj 



T'kjfik ) 



(A35) 



hok.k = \ ^hkkfl this becomes 



Making use of the gauge condition Eq. (All I, written as 
) this 

IGa- 



V^hoj — - —^hkkfij + hkkfij — hkj.ok 



We note that hkk = GGAe^^"^"?/ and h, 
IhkkMj, obtaining 



= -8nGe-^'l">pv^. (A36) 

kkfij — hkjflk = 



-i(v%,+6G.(|-lf)e-*.f,,, 

= -SttGc-^'^Vw^. (A37) 
Using the definition of the potential Vj and the relation 



Eq. ( A2 ) , we obtain 



= 2Ge"S'^"VVj. (A38) 



Finally, with the aid of Eq. ( A3 1 , we rearrange and solve 
for /loj 



'2_ IGa^ 
3 2 G 

(A39) 

Now to the spatial part of the vector. Eq. (12 1 to 
0(1.5) is 

K (^^A^ - + V'ho, - ho,,,, + hl%^ 

= -8TTGe~^'''°pv\l - e-^^«) (A40) 

Taking the covariant divergence of the vector equation, 
we obtain to 0(1.5) 



K . 



(1) 



2 "00, Oil 



-87rGe-2'^''(pw'),,(l-e^^'^«) (A41) 



using relations ( A24|A1 1 this can be written as 



KG 



-24>„y2^^^ = 2Ge'^*"'V^F,^,(l - 6""*'°), (A42) 



-200 rz^ 



400 \ 



or, with Eq. (A4| and canceling some terms on both 
sides, 



K 



-(7,o = -2(l-e-4*«)f/,o 



l-K/2 

From the requirement that U,o 7^ 0, we obtain 



l-X/2 



-2(1-6-'"^°) 



(A43) 



(A44) 



which links the cosmological value of the scalar field to 
the coupling constant of the vector field as follows 



' 4 \2^K 



(A45) 
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If the constraint Eq.(37l holds, it is possible to solve 
for the vector field to 0(L5) and obtain the PPN metric 
and the preferred frame coefficients, and this we proceed 
to demonstrate. Since the covariant divergence of the 
vector equation does not constrain we are free to 
set its value as we please. We can regard it as a C/(l) 
gauge transformation, with A" — > A" + (p'°'. Although 
the action is not fully U{1) gauge invariant because of the 
lagrange multiplier term, to 0(1.5), the relevant order for 
the preferred frame PPN parameters, the vector equation 
is gauge invariant. 

We thus set, for convenience 



With the definition of the potential Vi and 



we get 



K[A'- Ae-^'t'°GV, + -e-'^'>"'GK{Vi - W^) 



2e-20oG.i/^(l_g-4^.o) (A50) 



simplifying the vector equation 

^/iQi — ht^jji + 
= -STrGe-^-^Vw^l 



K ( V^A' 



2 ""00,20 



(A46) 



(A47) 



The spatial part of the vector field is then given by 



A' = e-^*" ( 2G ( 2e''^'^" + 



or, with Eq.(Al7l 



(A51) 



from Eq. (|A35p we have 



Oi 



yielding 

K ( V^A'' + 167rGe" 



= -8TTGpv'e~^'^°{l 



2 "-00,10 
-400 ^ 



(A48) 



Since h, 



(1) 

00 



■^hjj, we have 



K ( V^A* 



167rGe-^'^V«' " 
-STrGe-^-^Vw^l 



2 ^jj.Oi 



(A49) 



A' = 2Ge-2^° f f 2e-4^° + | 



{V, -Wi)\. (A52) 



This result satisfies A'j = if relation ( 37 ) holds 



We now combine Eqs. ( A39 1, ( A52 1 and ( A9 1 to obtain 



go^ = 2G(^^^^{V,-W^)~2e-*'^'>V^+ (A53) 



(1 - e-4^°) ( 



40o^ ( 1+{^K-I)e""f« 3 



K 



This remains to be converted to asymptotically 
Minkowskian coordinates using prescription ( A21[ ), with 
the potentials rescaled accordingly Vi = e^^^°Vi and 
W, = e-'^'I'oWi, yielding 



50. = 2G ( ''(2^!^)fV . - W,) - 2% + (e^^o - l) 



1 + (2ii: - l)e-'*^« 



K 



2-K 



(A54) 



In the GR limit, when K = k = (t>f) = Q (the limit 
should be taken before the division by K in the solution to 
the vector equation), and in units such that Gat = G = 1, 
this reduces to 



7 



1, 



(A55) 



as expected. 

Here we see again, that it is not possible to take K = 
without having 0o = Oi since then the physical metric 



would diverge. 



5. to 0(2) 



We have now arrived at the most complicated step- 
calculation of goo to 0(2). To this end, we should solve 

(2) 

for hf^f{ from the temporal part of the Einstein equations 



Eqs. ( A5 I, and for (j)^ ' from the scalar equation Eq. (Ill 



to 0(2). We start with /i, 



(2) 



12 



K 



,(1)|2 



,(1) 



-^^^'io - ^ Un,oo -2(1-^] h,o^,o 



(A56) 



All the terms in the second line disappear by virtue of 
the gauge conditions Eqs. (Alll,(A46l, leaving us with 



K 



2^(2) 



00 



-AnGe 



p I n + 3- + 2e-*'^«i;2 
P 



(A57) 



The left hand side of the equation yields, similarly to Using Eqs. (A31 1, (A25l we see that 
GR 



G 



2G 



K 



(v^^o^ + 2 {Gk? e-4'^"V2[/ - 8 (G^)' e-^^^V^^j 



2A;Ge"2'^V0^^^ = 2A:Ge"*"*° ( G^ + ^ ) C/p, 



fcG 



For the right hand side we use the definitions of the po- 
tentials $i,$2, ^"3, ^"4, together with Eqs. ( |A24[ ) , ( |A25p and with Eq.( |A23| this is 
to obtain 



(A60) 



-2e-2^« ( Gk 



4^ ) *2 + ^'3 + 3$4 



All in all, /iqq is given by 



/ijfc^Jj) - 2fcGe^20Op^(i) ^ 2kGe-^'''°GNUp. (A61) 



We now easily recognize the potentials <i>i,<i>2, <&3, <&4 and 
can write 



,(2) 
'00 



2e-*^« {Gxf + 4GKe-'**«$i+ 
8 (Ga-)' - 4Ga (gk - ^G^ ^ e-4*°<i>2 
2GKe-2'^''$3 + 6GAe-2'^°$4 



which simplifies to 

4? =2GAe~^"*° (-G/^e-^-^oC/^ + 2e-''"*°$i 

+2e-2'^°GAr«'2 + *3 + 3«'4) (A58) 

Lastly, the scalar equation to 0(2) is 

V20(2) - /ij-fe0(i) + 2fcGe-2^V0(i) = 
A;Ge-2'^° (2pe-^'^" ^_ + 3p) (A59) 



0(2) 



fcGe-^-Jio 
47r 



(2e-'''^«$i + 2e-2'^''GAr$2 + ^'3 + 3*4) 

(A62) 



We now have all that is required to find 500 to 0(2). 
We first write the contribution of the 0(2) term 



hg^ ~ 20(2) - 2(0(1) )2 + 20(1)41) = 2GAre-20o 

(2e~^'^«$i + 2e-2'^"GAr$2 + *3 + 3$4 - e''^'*'" G nU"^) 

(A63) 



Combining all orders of goo i we get 
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goo 



=200 



$3 + 3*4)) 



(A64) 



Transforming coordinates according to ( |A21 1, rescaling 
the potentials to the new coordinates and choosing units 
such that Gat = 1, we finally obtain 



and 



5oo 



-1 + 2U- + 4$i + 4$2 + + 6^4, (A65) 



exactly as in GR. From here we can instantly read that 
03 = Ci = C2 = C3 = C4 — 0; expected for equations 
of motion derived from a lagrangian. We also see that 
/3 = 1 and ^ = 0. 



APPENDIX B: PPN PARAMETERS FOR TEVES 
WITH jETHER type ACTION 

All the setup for the calculation of the PPN parame- 
ters remains as in Section |A] except for the gauge condi- 
tion for the time-space component of the Einstein metric, 
which we shall recalculate. ^* and are as in Section 
A 1 We follow the same steps as in the calculation for 



regular TeVeS , and for brevity we only show the relevant 
equations and altered results. The main difference will 
be in the calculation of goi, since the divergence of the 
vector equation will now constrain the spatial divergence 
of the vector field, and not the coupling parameters of 
the theory (see Ref. [50], Section 5.4 for details). This 
will enable us to calculate goi unambiguously, and obtain 
the preferred frame parameters. 



(B5) 



with Gy replacing Gk- The physical metric in asymp- 
totically Minkowski coordinates is 

goo = -1 + 2GnU , gij = Sij ,got = (B6) 

with Newton's constant redefined as 



Gtv = G 1 - 



K + K+- lU 



k 

An 



(B7) 



2. g^j to 0(1) 



The result here is again the same as in Section |A 3[ 
with Gu replacing Gk 

hij =2Gye-'^'t'"U5ij, (B8) 

and in asymptotically Minkowski coordinates, with 
Gjv = 1, 

9^] = (1 + 2U) . (B9) 

Comparing this with Eq. ([T6| we see that for TeVeS with 
iEther type action, the PPN parameter 7 is still unity. 
This is a promising result for the modified theory, since 
solar system data [H] strongly constrains the value of 7 
to be near unity. 



1. 500 to Oil) 
The temporal component of the Einstein equations 



(A5) with the altered stress tensor Eq.(40l, is to 0(1) 



2'^ 



K + K^- Ka 



2^(1) 



00 



47re-2'^°Gp (Bl) 



yielding for h. 



(1) 



feW^2e'^^"G(l- ^ + ^+'^^ 



3^' 



For convenience, we define 



G„ EE G 1 



K + K+-Ki 



(B2) 



(B3) 



In fact, this is the quantity that replaces Gk in the cal- 
culation. With the scalar equation unaltered, we have, 
as in |A2l 



'■00 



2G^,f7 



(B4) 



3. goj to 0(1.5) 



The Lagrange multiplier A to 0(1) is now given by 
1 
2 



A = ^ (X - K+) V^^o - 167rGpsinh20o (BIO) 



The time-space component of Eq. ( A5 1 with the new 



stress tensor Eq. (40 1 to 0(1.5) is then 



- (V^/ioj — hak,kj + hkk,oj — hkj.ok) ^hkk.oj' 

^ {K2 + 2K+) - {K2 + K+) A% - K+A^,, = 



and the vector equation to 0(1.5) is 



(Bll) 



{K + K+- Ki) /i^o^o ■ + K {hoj^kk - hok^kj) + 



1 
2 

^K2hkkfij + K+hjk.ko + {K+ + K2-K) A%+ 

+ {K + K+) ^ -87rGe-2-#'o(l _ e-*'^«)pwJ' 

(B12) 
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These two equations should be solved to yield hoj and which simplifies to 



. We still need two ingredients to be able to solve 
the equations A-'j, which will be determined from the di- 
vergence of the vector equation, and hok.k, whose value 
should be set from the demand that 500 be in the stan- 
dard PPN gauge. In the meanwhile, we write 



l-Qk.k 



-200 



a i3 + Q)U,o, 



(B13) 



A'. 



2K^ 

3K: 



(2(1 



-K2 
K + 3K2 



2-{K + K+~- Ki) 



(B16) 



and we will solve for Q in the next step of the calcula- 
tion. This is inspired by the standard GR gauge condi- 
tion hok.k = 3?7,o, since in the GR limit, Gy ^ G = 1 
and = 0. The divergence of the vector equation is 



and finally, with Vjj — —U,o = -^Xfljji we have for the 
spatial divergence of the vector field 



^{K + K+- K4) 4o,ow + + K+hjk,k]0+ 



+ {K+ + K2)A%^+K+A^,,^ 



(B14) 



Using results from previous steps, this can be written as 

-i(3i^+ + K + 3K2- K4)e-^^'>GyX,0nkk 
+ {2K+ + i^2)^',-fefe = 2Ge-2^"(l - e-^*")Vj,okk (B15) 



■ _ Ge-^'t'" 



1 



ZK+ + K + 3K2- K4 



2-{K + K+-Ki) J 
(B17) 

The addition of extra kinetic terms in the vector action 
has removed the U{\) gauge invariance, and enabled us 
to determine unambiguously the spatial divergence of the 
vector field. 



X, 



We may now substitute the vector divergence Eq.(B17l 
and the metric gauge Eq.(B13l into Eqs. ( B12|B11 1 and 
solve for and /loi- Doing so, we obtain 



4Ge-2'^°(Xe-4'^o + K+)Vi 
2KK+ -~{K + K+) 

Ge-^'f'« {2K{2 -{K + K+- X4))e-^^» + {2KK+ - {K + K+))Q + K+-'iK + 2K+{Ki -K+ + 2K)) x,o^ 

{2-{K + K+- Ki)){2KK+ - {K + K+)) 

(B18) 

2Ge-^^"((l-2X)e-4^»-l) / {K + {2K ~ 1){K+ + K2))e-^^» 

2KK+-{K + K+) * ' \{K2 + 2K+){2KK+^{K + K+)) 
4(1 ~ K)Kl - 2{Ki + K2{3K - 2) + 1)K+ + 2(1 + 2K2)K - [K^ + 2)K2 \ 

{K2 + 2K+){2KK+-{K + K+)){2-{K + K+~K4)) ) ^ ^> 



We still need to find Q; to this end, we proceed to the 
last step in the calculation. 



4. 500 to 0(2) 



(2) 

The equation for h^Q is 



15 



]^{K + 3^2 + iK+ - Ki) -0 - i (^1 + + /j,,,oo + ^ (2 - (X + X+ - K^)) hjo,jO 
: -47re-2'^>oGp (/i^ + 20(i)) + 47re-2'^''Gp (^H + 3^ + 2e-^'l">v^^ (B20) 



I 

All terms on the second line are proportional to x.oo; requirement 
therefore the standard PPN gauge is obtained from the 



{K + 3K2 + 3K+ -Ki)A^.--il + -K2 + K+ h,,^o +-{2-{K + K+- K^)) /i^oj = (B21) 



or 



Ge-^-^o {3K+ + K + 3K2- K4 
2K+ + K2 



1 - e 



-400 



3K+ + K + 3K2- K4, 
2-{K + K+- K4) 



(B22) 



Combining this with (B13l, we obtain for Q 



?>K+ + K + m2-iu . 

^ 2K^ + K, 



2-4X4 



2-{K+ + K- K^) 
(B23) 



Using this value of Q, in quasi-Cartesian coordinates, 
with the potentials appropriately rescaled, the time-space 
component of the physical metric is given by 



9oi 



'2 ((2 + 3K2 + 2if+)e2'^n -{2-{K + K+- K^)) e-^^o) gi^h (20o) 
^ {K2 + 2K+){2-{K + K+-Ki)) 

2 ((e^-^o - (1 - 2X)e-2'/>o) sinh (20o) - (i^+e*"^" + K)) 
2KK+ -~{K + K+) 

{{Kj -K + 3K+){2 + 3K2 + 2K+)e^*« + {m+ + 3K2 + K- Ka){2 ^ {K + K+ - Kj))) ' 

{K2 + 2K+) {2-{K + K+ - Ki)f _ 
' {{K4 -K + 3K+){2 + 3X2 + 2K+)e^'t'° + {3K+ + 3K2 + K- Ki){2 - {K + K+ - K4))) 
, {K2 + 2K+) {2-{K + K+- Ki)f 

2 ((2 + 3K2 + 2K+)e^'l>« - (2 - (X + K+ - K^)) e-20o) ginh (2^o) 

{K2 + 2K+){2-{K + K+-Ki)) 
2 ((e^-^o - (1 - 2if)e-2<^«) sinh (2(/)o) - (i^T+e^^" + X)) ' 



2KK+ -{K + K+) 



GWi (B24) 



In the standard PPN gauge, for TeVeS with iEther 

(2) 

type vector kinetic term, the equation for /iqq is the 
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same as for simple TeVeS, Eq.(A56l, except for the re- 
placement of Gk by Gy. Since the scalar equation 
remains unchanged, we can immediately write the re- 
sult for goo 5 after coordinate redefinition and appro- 
priate rescaling of the potentials, and this time with 



Gjv-G( (1- 



47r 



1 



once again we have 013 = = ^2 = Ca = C4 = Oi ^ ex- 
pected for equations of motion derived from a lagrangian, 
and also f3 = I an d — 0, as for simple TeVeS. We can 
now compare Eq.(B24l to the standard form Eq. (151, 



and extract the preferred frame parameters 



5oo 



-1 + 2U- 2U^ + 4$i -I- 4$2 + 2$3 + 6$4, (B25) 



^ G ((e^^o - (1 - 2j^)e-2^°) sinh (20o) - jK+e^^" + K)) \ 

2KK+-{K + K+) J ^ ' 

^ oi AG ((2 + 3K2 + 2K+)e^'i>» - {2 ~ {K + K+ - K4)) e-^-f") sinh (20o) 

2 {K2 + 2K+){2~{K + K+~Ki)) ^ ^' 

2G {{Ki -K + 3K+){2 + 3K2 + 2K+)e^^° + {3K+ +3K2 + K- K^){2 - [K + K+ - K^))) ^ , 

^ h o (rJZo) 

{K2 + 2K+){2-{K + K+-Ki)f 



In the above, 

47r(2 -{K + K+) + K4) 



in + k{2 -{K + K+) + Ki) ■ 
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